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We investigate the diffusion limited aggregation of particles executing persistent 
random walks. The scaling properties of both random walks and large aggregates 
are presented. The aggregates exhibit a crossover between ballistic and diffusion 
limited aggregation models. A non-trivial scaling relation ^ ~ ^^-^^ between the 
characteristic size ^, in which the cluster undergoes a morphological transition, and 
the persistence length i, between ballistic and diffusive regimes of the random walk, 



'^ is observed. 
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lO 1 Introduction 
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,__! Aggregation processes are outstanding phenomena that have been subject 

L| of theoretical, experimental and simulation investigations in many fields of 

•^ knowledge [T1I2] . The diffusion limited aggregation (DLA) model proposed by 

*^ Witten and Sander [3] is the most studied theoretical aggregation process. In 

^ this model, free particles are released, one at a time, far from a growing cluster 

and perform successive random jumps while they are not adjacent to the clus- 
ter. If a free particle encounters any particle of the cluster, it irreversibly sticks 
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in the corresponding position and becomes part of the cluster. Notwithstand- 
ing the rule simplicity, the DLA model produces clusters exhibiting complex 
scaling properties P|I5|5|7] . Another simple and well studied growth process 
is the so called ballistic aggregation (BA), in which the free particles move 
ballistically at randomly chosen directions and obey the same sticking rules as 
the DLA model. The BA model produces asymptotically homogeneous radial 
clusters [8] whose active (growing) zone is described by Kardar-Parisi- Zhang 
(KPZ) universality class [9J, as observed in the planar version where the flux 
of particles is normal to an initially fiat substrate p^fTT] . 

Several generalizations of the DLA model have been studied [1] with special 
attention in those where the free particles perform random walks with a drift 
[T^13fl4IIT^ . In these models, the clusters undergo a crossover from a DLA to a 
BA scaling regime as the number of particles increases due to the trajectories 
become essentially ballistic at asymptotic large scales. In another group of 
models, in which the free particles perform long random steps of a fixed length, 
an inverted transition from BA to DLA is observed p^JITT] . In this case, the 
trajectories are essentially random at asymptotic large scales, independently 
of the step length. A simple characterization of the morphological transition is 
done by relating the mass M of an aggregate with its radius r. This quantity 
has two asymptotic scaling regimes, M r^ r"^ and M ~ r^'^^, related to either 
the DLA or BA models, respectively. 

Aggregation models with drift or long steps can be used to investigate systems 
with a long mean free path as, for example, the aggregation of methane or 
ammonia molecules in superfiuid helium [18] . In particular, an off-lattice DLA 
model with long steps of length 5 was recently applied to the aggregation of 
molecules in superfiuid media [T^. In addition to this example, biased diffusion 
and super-diffusivity have also been reported in several systems ranging from 
adatom surface diffusion p^f20] to cell migration [2T1I22] . The relevance of 
modelling the aggregation with this kind of random walk is strengthened due 
to biased diffusion may be misleadingly confused with super-diffusivity [23] . 

The representation of a large mean free path by a ballistic movement is a 
simplifying hypothesis since weak perturbations can promote slight deviations 
of the ballistic path resulting a biased diffusion at short scales. Motivated by 
the wide availability of systems for potential applications, we investigate a 
two-dimensional aggregation process where the particle trajectories consist 
of random walks, with the direction of a new step limited by an angle 80 
in relation to the previous step direction. This correlated walk, which was 
formerly investigated in other contexts [M|25l)26] . has a drift at short scales 
and becomes a random walk at large scales. Also, an aggregation model using 
this random walk was formerly investigated in a small cluster size limit and a 
fractal dimension depending on the 5o parameter was reported [15j . Actually, 
this dependence is a finite size effect since this kind of trajectory lead to a BA 



to DLA crossover, as demonstrated in the present work. 

The paper is organized as follows. Model and simulation strategies are pre- 
sented in section [2j The trajectories are numerically investigated and an an- 
alytical approach for the trajectory crossover is developed in section [3j In 
section |4| the simulations of large clusters are presented and a scaling analysis 
of the morphological transition is developed. Finally, some conclusions and 
prospects are drawn in the section |5} 



2 Model 



We perform two-dimensional off-lattice simulations with an initial cluster con- 
sisting of a single particle of diameter a (a seed) stuck to the origin. Free par- 
ticles are sequentially released, one at a time, at a circle of radius r/ 3> r^ax 
centred at the seed, where Vmax is the largest distance from the seed of a par- 
ticle belonging to the cluster. The free particles follow a persistent random 
walk [24j, in which the position of the nth step is given by 



Xn =Xn-l +aCOS(j)n 

?/n = l/n-i + asin0„ (1) 

and the direction of the nth step depends on the preceding one as 

4>n = 4>n-l + Vn- (2) 

The white noise rjn is a random variable uniformly distributed in the interval 
{—6g/2, 60/2) and the parameter 6g limits the next move direction inside an 
angular opening of size 6g centred in the direction of the previous step. 

The trajectory becomes ballistic for 5^ — > and random for 6$ — )■ 2tt implying 
in the BA and DLA limit cases, respectively. The direction of the first step 
is chosen at random. The trajectory is stopped whenever the particle visits 
a position adjacent to the cluster where it irreversibly sticks. Finally, the 
particle is discarded whenever it crosses a distance r^^ ri = Vmax + A. The 
procedure is repeated up to the cluster reaches N particles. The previously 
introduced variables r; and r^ must be as large as possible, but computational 
limitations restrict their values. For the DLA simulations, A can be a few 
particle diameters [27] • However, for the BA limit, A cannot be small due to 
distortions caused by shadow instabihties [28|29] . The BA clusters are free 
from shadow effects when A > 600 and thus A = 800 and r^ = 10 x r^ were 
adopted in all simulations. 

Since the asymptotic behaviour be typically reached only for large clusters. 
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Fig. 1. (Colour on-line) Angular distributions of the first passage at a distance r 
for a persistent random walker starting at the origin. Without loss of generality, we 
set (pQ = 0. Symbols represent the stochastic simulations and solid lines non-linear 
regressions using Eq. Q. 

an efficient algorithm is required. We adopted the strategy where particles 
execute long steps with the proper angular distribution P{4>) if they are not 
close to the cluster [12] . A closed expression for the angular distribution of a 
long jump is not known a priori. So, we numerically computed the probability 
of the first passage in a circle of radius r occurring at an angle for a trajectory 
started at the origin. Figure [T] shows some angular distributions for different 
values of 6g and r obtained with up to 10^ independent trajectories. If r and 
6g are not too large, the distributions are very well fitted by 



P(0) = Aexp(5cos0) + a 



(3) 



Parameters A, B and C (actually only two of them are independent due to 
the normalization) were determined using non-linear regressions that are also 
shown in Fig. [1} Cluster growth optimizations were implemented using steps 
of size 16a in the large empty regions nearby the cluster and three step sizes 
(16a, 100a and 200a) were used far from the cluster. More details about long 
step and off-lattice optimizations are available elsewhere 



3 Trajectory scaling properties 



Figure [2] shows two stages of a trajectory corresponding to short and long 
times. Roughly, the trajectory consists in a drift around the initial direction 
until a large angular fluctuation be selected, as illustrated in the left panel of 
Fig. [2j The smaller angular opening the rarer the probability of a large flue- 





Fig. 2. A typical trajectory for 5e = 10° after 10^ (left) and 10^ (right) steps. 

tuation be selected. However, these rare events unavoidably occur for a suffi- 
ciently long time if ^g 7^ and the trajectories become erratic at large scales, 
as shown in the right panel of Fig. |2} In order to investigate this crossover, 10^ 
independent trajectories were simulated. In Fig. 3(a) the root-mean-square 
(RMS) displacement A is show for several values of 5g. One can clearly resolve 
a crossover: the RMS displacement scales as A ~ t at early stages and crosses 
over to a scaling A ~ t^/^, as expected in a transition from ballistic to diffusive 
trajectories. 

The characteristic crossover time depends on the angular opening as r ~ 
Sq"^. This result can be proved with the following reasoning. Without loss of 
generahty, we assume 0o = in Eq. ^ resulting 



0t = Zl ^fc' 
fc=i 

where 77 is a white noise with (r/) = 0, {rf) oc (5|, and (■ ■ ■ ) means statistical 
averages. Obviously, the turning angle is a random variable with (0j) = 0. 
Applying the central limit theorem [SU] , in which the sum of t identical and 
independent random variables with average {rj) = and variance (r^^) con- 
verges to a normal distribution with variance af = t{rf), the probability of 
a large turning angle be selected is exponentially negligible for ctj ^ 1 and 
appreciable for crj 3> 1. Assuming a crossover time given by 



(Tr 



T{r]^) 



we immediately find out the scaling law r ~ 6^'^. Finally, the characteristic 
RMS displacement at the crossover is proportional to r since the trajectory 
is still approximately ballistic at the crossover. These scaling properties A can 
be expressed in the following scaling ansatz 



A(t, Sg 



tail 



(4) 



where the scaling function g{x) has the properties g{x) ~ const, for x ^ 1 
and g{x) ~ x~^'^ for x ::^ 1. As one can see in Fig. 3(b) this ansatz collapses 
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Fig. 3. (Colour on-line) (a) RMS displacements for distinct angular openings. 
Dashed lines represent scaling laws A ~ t and A ~ t^' ^. (b) Collapse of curves shown 
in Fig. |3(a) using the scaling ansatz given by Eq. Q. The dashed line represents 
the asymptotic behaviour of the scaling function g{x) ~ x^^'"^ . 



the data shown in Fig. 3(a) onto a universal curve if we plot A/t against t x (5g 



demonstrating the ansatz correctness. 

Before analysing the aggregates obtained using the persistent random walk, we 
shortly draw a comparison with the results reported by Tojo and Argyrakis for 
this persistent random walk [23], in which a scaling r ~ ^~q^'^^ was numerically 
obtained in stochastic simulations. This exponent, that is slightly different 
from ours, results from the crossover times estimated using the intersections 
between early and long time scaling laws of A versus t. We instead deduced the 
scaling behaviours and the theory perfectly matches the numerical simulations 



as shown in Fig. 3(b) 



4 Aggregate simulations and scaling analysis 



This section is devoted to the simulations of a DLA model using the persistent 
random walk defined by Eqs. ^\ and (|2]). Figure 111 shows different stages of 
two clusters obtained for angular openings ^q = 40° and 60°. In both cases, the 
change of the cluster morphology is evident. At early stages, the aggregates 
have dense branched morphologies (low density of inner voids) that resemble 
the BA patterns while sparse branches, the hallmark of DLA clusters, are ob- 
tained after a characteristic number of particles. This behaviour qualitatively 
confirms the morphological crossover predicted with the trajectory properties: 
the smaller Sq the longer the crossover time. Figure |5] shows clusters of the 
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(d) A^ = 10^ (e) N = 10^ (f) N = lO'^ 

Fig. 4. Different growth stages of clusters generated for 6g = 40° in (a)-(c) and 60° 
in (d)-(f). The number of particles are show below each panel. 

same size obtained with distinct 6g values. The simulations were stopped when 
the clusters reached a circle of radius r = 1000a. A transition from dense to 
ramified structures is observed and a fractal dimension dependent on Sg can ef- 
fectively be measured in agreement with Ref . p^ . However, it is important to 
notice that the dense morphologies are transient and the patterns unavoidably 
become ramified with the DLA fractal dimension {D ~ 1.71) for asymptotic 
large clusters. 

The quantitative characterization of the crossover was done using the mass- 
radius method [8]. Clusters were grown until they cross a circle of radius 
5 X lO^a. In order to perform statistical averages, 100 independent samples 
were simulated for each investigated 6g value. The mass-radius method consists 
in determining the number of particles inside a circle of radius r centred in 
the seed. In general, the mass scales as M{r) ~ r^, where D is the fractal 
dimension. We have D ~ 1.71 and D ^ 2 for asymptotically large DLA and 
BA clusters, respectively [81IT2] . 



Figure |6] shows the mass-radius curves corresponding to different values of 6g. 
The curves have a crossover between the scaling regimes M ~ r^ at small 

at large radius. The mass is described by a scaling form 



radius and M 



,1.71 



where /(x) 



X 



M{r,Se) = 
for X <^ 1 and /(x) 



-22 



/ i^er) 



(5) 



X 



1.71 



for X ^ 1. This scaling ansatz 



supposes that the crossover radius scales as ^ ~ dg"". It also assumes a con- 
straint M(^) ~ ^^ due to the clusters still are dense close to the crossover, 
explaining the exponent 2z in Eq. ([s]). Determining the crossover radius as 
the intersection between the scaling laws Mi ~ r^ and M2 ~ r^-^^, we ob- 
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(b) 5e = 20° 




(c) 5e = m° 



(d) 5e = 120° 



Fig. 5. Clusters of a fixed size r ^ 1000a obtained for distinct Sg values. The number 
of particles varies from 10^ for the most ramified to 10^ for the densest cluster. 
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Fig. 6. (Color on-line) Mass-radius curves for three 5g values. Dashed lines are power 
laws with exponents 2 (smaller r) and 1.71 (larger r) included as guides to the eyes. 
The inset shows the mass M* and the characteristic radius ^ at the crossover as 
functions of Sg. Power law regressions with exponents —5 and —2.5 are also shown 
in the inset. 

tained a scaling law very close to ^ ~ ^g^'^- We also evaluated the mass at the 
crossover and the relation M(^) ~ ^"^ was confirmed. The crossover mass and 
radius against 6g are shown in the inset of Fig. [6] 



In Fig. [Tithe mass-radius curves are collapsed using the scaling ansatz given 
by Eq. pi) with the exponent z = 2.5. As can be seen, an excellent collapse 




Fig. 7. (Colour on-line) Collapses of the mass-radius curves for 6g = 10°, 20°, 40°, 
60°, 120°, and 240° using the scaling ansatz given by Eq. ([s]). The solid lines are 
scaling laws r^ and r^'^^ corresponding to the asymptotic behaviour of the scaling 
function f{x). 

onto the scaling function f{x) is obtained. For the general case with a drift 
quantified by a parameter b, the mass against radius can be written as [TTj 



M{r; h) = h'^f 



(6) 



wliere the scaling laws ^ ~ 6^ and M(^) ~ 6" are assumed for the crossover 
radius and the corresponding mass, respectively. The scaling function f{x) 
yields the asymptotic regimes of M{r): For a DLA to BA transition, we have 
f{x) ~ x^'"^^ for X ^ 1 and f{x) ~ x^ for x ^ 1 while the scaling exponents 
are exchanged to f{x) ~ x^ for x ^ 1 and /(x) ~ x^'^^ for x ^ 1 for a 
BA to DLA transition. We can easily relate the scaling hypothesis rtsl) with 
the scaling ansatz given by Eq. (|6| for the aggregation of particles performing 
long steps of size £ = 6 . It was shown in Sec. |3]that the characteristic length 
in a persistent random walk scales as £ ~ 5^^. Replacing 5q ~ 1/yJ in Eq. 
(|5| and comparing with Eq. ([6]), we find the exponents a = z = 2.5 and 
C, = z/2 = 1.25. These exponents are consistent with those obtained for long 
steps trajectories in Ref. [17] indicating a universality class. 



5 Conclusions 



In this work we investigated the scaling properties of a diffusion limited aggre- 
gation (DLA) model where the particles follow locally persistent trajectories. 
The bias is introduced by limiting the direction of a new step into an angular 
opening Sq in relation to the preceding step direction. We show analytically 



and numerically that this trajectory is effectively ballistic below a length scale 
i r^ 6g'^ and becomes random at large scales. This result improves the nu- 
merical estimate of the scaling law reported formerly by To jo and Argyrakis 

m- 

Large scale simulations show that the aggregates undergo a morphological 
transition from a ballistic aggregation (BA) regime at small sizes to a DLA- 
like ramified morphology at asymptotically large scales. It is important to 
mention that a previous analysis of this aggregation model [15], which was 
done in a small cluster size limit, states a non-universal fractal dimension. 
However, we have shown that the asymptotic behaviour of the clusters is 
described by the DLA fractal dimension. Finally, the crossover between the 
BA and DLA scaling regimes is related with the angular opening by ^ ~ ^7^^^ 
indicating a non-trivial relation between the characteristic size of the cluster 
in the crossover and the persistence length of the random walk. 

The persistent random walk has interesting features that can be associated 
to a number of physical systems [18,22,26,31,32j . A small angular opening 
may represent a diffusion in a medium where the scattering is weak except 
by rare but no negligible fluctuations. An important example is the wander- 
ing/aggregation of molecules in superfiuid helium [18] where the particles may 
have a mean free path ranging from a molecule diameter to micrometers de- 
pending on the proximity to the critical temperature. This kind of diffusion is 
also observed in cell migration pTP2] , a chemotactic process in which the cells 
diffuse following a gradient-mediated chemical signalling. A relevant property 
of the trajectory is a characteristic mean free path or persistence length that 
emerges naturally as we have demonstrated with the central limit theorem |3Dj . 
In the cited examples, the aggregation and the consequent formation of clus- 
ters were subjects of recent experimental interest [TS|33ll34j . Consequently, the 
present theoretical analysis can be extended to future applied investigations. 
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